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1.  Introduction 


Mercury  cadmium  telluride  (HgCdTe)  alloy  is  of  great  importance  in  sensing  radiation  from  the 
near  infrared  (A,c  ~  1  pm)  to  the  very  long  wavelength  infrared  (Xc  -  15  pm).  Much  of  the 
HgCdTe-related  research  and  development  work  is  carried  out  for  cooled  operation.  Intrinsic 
carriers  play  a  dominant  role,  especially  at  long-wavelength  (8  pm  to  12  pm  cut-off)  material 
near  ambient  temperatures  due  to  high  thermal  generation  of  carriers.  This  results  in  low 
minority  carrier  lifetimes  due  to  high  Auger  recombination  processes.  Consequently,  this  low 
lifetime  at  high  temperatures  results  in  high  dark  currents  and  high  noise.  Cooling  is  one  means 
of  reducing  this  type  of  detector  noise.  The  challenge  is  to  design  photon  detectors  to  achieve 
background-limited  performance  (BLIP)  at  the  highest  possible  operating  temperature,  with  the 
greatest  desire  being  operation  close  to  ambient  temperature.  Condition  of  BLIP  is  defined  as  the 
noise  being  entirely  due  to  the  external  photon  flux  and  not  due  to  the  detector.  Hence,  BLIP 
condition  depends  on  the  incident  photon  flux,  temperature,  detector  dark  current,  etc. 

The  trick  of  the  trade  is  to  get  high-operating-temperature  (HOT)  detectors,  which  are  being 
designed  to  operate  at  -AT  K  higher  temperature  than  existing  detectors,  but  with  performance 
matching  at  low  temperatures.  Another  option  is  to  operate  these  detectors  at  the  same 
temperature  as  the  existing  detectors,  but  with  higher  performance. 

A  higher  operating  temperature  would  result  in  several  advantages  to  an  infrared  imaging 
system — mainly  a  reduction  in  power  and  weight  requirements,  thus  lowering  mission  costs  and 
allowing  infrared  systems  to  be  deployed  on  various  new  platforms.  This  work  would  have 
significant  impacts  on  several  mission  areas  including  the  Missile  Defense  Agency,  the  Army, 
and  other  platforms.  This  technology  will  increase  IR  detection  range  and  spatial  resolution  of 
imagery,  and  reduce  the  carry-on  weight  of  an  individual  Soldier,  and  has  a  significant  impact  on 
the  size,  weight,  and  power  (SWaP)  of  the  overall  IR  imaging  system. 

Limiting  condition  is  the  existence  of  high  intrinsic  carriers  near  ambient  temperatures. 
Calculated  dark  current  density  of  a  10  micron  cut-off  (at  78  K)  detector  versus  temperature  is 
shown  in  figure  1  together  with  noise  currents  calculated  for  different  photon  fluxes.  As  can  be 
seen  at  80  K,  dark  current  density  is  in  the  range  of  high  106  Amps  per  cm2  and  exponentially 
increases  as  the  temperature  increases.  At  150  K,  dark  current  density  is  in  the  low 
10  “  Amps/cm"  range.  The  choice  of  detector  architecture  determines  the  type  of  dark  current 
that  impacts  detector  performance.  Under  optimized  conditions,  dark  saturation  current  density 
of  a  state-of-the-art  photovoltaic  (PV)  detector  can  be  written  as  (see  later  sections  for  the 
derivation), 
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where  Jodiff  denotes  diffusion  current  density  in  the  absorber  layer,  rp  denotes  hole  minority 
carrier  lifetime,  Nd  denotes  n-type  donor  concentration  level,  rij  denotes  the  intrinsic  carrier 
concentration,  and  d  is  the  thickness  of  the  absorber  layer. 

As  shown  in  figure  lb,  the  detector  is  not  at  BLIP  at  100  K  for  photon  flux  of  less  than 
1015  ph/sec/cm2  but  shows  BLIP  greater  then  photon  flux  of  1015.  Hence,  dark  current 
suppression  is  the  key  to  increasing  the  operating  temperature  or  increasing  the  performance  at  a 
particular  temperature.  As  can  be  seen  from  the  previous  equation,  to  reduce  the  dark  current,  the 
thickness  of  the  absorber  layer  and  the  density  of  the  minority  carriers  need  to  be  reduced,  and/or 
minority  carrier  lifetime  has  to  be  maximized. 


Figure  1.  (a)  Calculated  dark  current  density  versus  temperature  of  a  HgCdTe  photovoltaic  device.  Cut-off  is 

10  microns  at  78  K  and  (b)  calculated  dark  noise  and  noise  due  to  photon  fluxes  versus  temperature.  Noise 
is  in  Amps  per  cm  per  root  Hz. 

HOT  technology  has  been  untenable  because  of  several  severe  difficulties  in  reducing  the  dark 
current.  The  fundamental  approach  that  was  used  to  increase  the  operating  temperature  was  to 
reduce  the  volume  of  the  diffusion  region  of  the  device  while  maintaining  the  performance  at 
low  temperatures.  Reducing  the  diffusion  volume  of  a  detector  can  be  achieved  by: 

a)  Reducing  the  active  layer  thickness.  This  approach  significantly  impacts  the  absorption  of  the 
photons  in  the  long  wavelength  spectrum  of  the  spectral  band  pass.  When  the  layer  thickness  is 
much  less  than  the  cut-off  wavelength  (dark  current  is  linearly  proportional  to  the  absorbing 
layer  thickness),  quantum  efficiency  (QE)  at  longer  wavelengths  will  suffer  drastically.  Figure  2 


shows  an  example  of  such  a  case  (dotted  curve).  To  overcome  this  problem,  one  has  to  have  a 
double  pass  (or  multiple  passes),  which  is  a  challenge  in  the  device  fabrication  process. 


Figure  2.  (a)  Dependence  of  absorber  thickness  on  quantum  efficiency  versus  wavelength  of  a  long  wave- 

HgCdTe  device  (Cut-off  0.c)  is  10  pm)  at  78  K.  (b.  Calculated  dark  current  density  versus  absorber 
layer  thickness  at  78K. 

b)  Reducing  the  junction  area  while  maintaining  the  optical  area  as  the  original  pixel  by  using 
micro-lens  technology  on  a  pixel  level  (3)  is  shown  in  figure  3.  Small  junction  area  will  increase 
the  device  impedance  and,  hence,  reduce  the  dark  current,  thereby  reducing  the  detector  Johnson 
noise  and  permitting  longer  integration  time  for  improved  signal-to-noise.  Additionally,  this 
architecture  will  reduce  the  influence  of  defects  intersecting  the  junction,  hence  increasing  in 
operability  (4).  One  of  the  challenges  is  difficulty  in  producing  micro-lenses  in  a  small  pitch 
format  on  the  backside  of  a  substrate. 
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Figure  3.  Cross-sectional  drawing  of  micro-lensed  hybrid  focal  plane  array  of  HgCdTe 
on  CdZnTe  substrate. 

c)  Recently,  DARPA  was  soliciting  an  innovative  device  architecture  inspired  by  nature  (5).  This 
is  based  on  a  moth-eye  nano-pillar-like  texture  to  reduce  reflection.  Additionally,  these  pillar-like 
structures  will  decrease  the  diffusion  volume  of  a  detector.  Hence,  photon-trapping  device 
architectures  give  the  same  sensitivity  as  a  conventional  device,  but  a  much  reduced  dark  current 
permitting  operation  in  the  higher  temperatures.  The  aim  of  the  DARPA  program  is  to 
demonstrate  these  nano-pillar  structures  as  an  inherent  anti-reflection  structure  (similar  to  moth- 
eye)  while  performing  its  function  as  a  detector,  possibly  for  HOT  applications.  One  of  the  issues 
related  to  this  architecture  is  the  hardness  of  surface  passivation  due  to  high  surface  to  perimeter 
area  (see  figure  4). 


Figure  4.  Schematic  example  of  a  one  pixel  pillar  like  structure 
with  16  sub-element. 

As  explained  earlier,  besides  reducing  the  diffusion  volume,  other  methods  of  reducing  the  dark 
currents  are: 
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1)  Increasing  the  minority  carrier  lifetime,  and 

2)  Reducing  the  minority  carrier  concentration. 

This  report  will  henceforth  describe  these  two  methods  of  reducing  dark  currents  of  a  given  PV 
device. 

The  dominant  intrinsic  recombination  mechanism  in  long-wavelength  HgCdTe  is  Auger 
recombination.  The  electron-initiated  Auger- 1  process  dominates  in  n-type  material,  whereas  the 
hole-initiated  Auger-7  process  dominates  in  p-type  material.  Intrinsic  Auger-7  lifetimes  are 
significantly  longer  than  intrinsic  Auger- 1  lifetimes  ( 1 ,  2).  Hence,  p-type  HgCdTe  absorber 
material  is  preferred  for  higher  operating  temperature  application  devices  (2). 


2.  Nonequilibrium  Mode  of  Operation-Auger  Suppression 


Recently,  Ashley  and  Elliott  (6-8)  proposed  a  new  class  of  infrared  sensing  devices  that  used  the 
phenomena  of  exclusion  and  extraction  of  thermally  generated  carriers  from  the  active  region 
under  reverse  bias,  thereby  suppressing  Auger  processes.  I  will  be  discussing  device  architecture, 
which  relates  to  exclusion  and  extraction  junction  in  HgCdTe  P+/p/N+  configuration.  Similar 
work  was  done  by  Elliot  (6)  in  indium  antimony  (InSb)  material  system.  In  this  report,  I  am 
discussing  a  P7ti/N+  HgCdTe  structure,  where  the  symbol  “ti”  indicates  a  drastically  reduced 
extrinsic  p-type  carrier  concentration  in  the  order  of  mid  10 15  cm-3;  P+  and  N+  denote  higher 
doping  density,  as  well  as  a  higher  energy  gap  than  the  photosensitive  base  p-region.  In  the  n- 
region,  the  minority  carriers  are  electrons.  When  the  minority  carriers  are  electrons,  a  gain  of 
approximately  60  times  in  the  minority  carrier  lifetime  is  obtained  (see  figure  8)  because  of  the 
lack  of  symmetry  between  the  valence  and  conduction  bands  in  LW-HgCdTe,  regardless  of 
doping. 

The  challenge  to  this  approach  is  to  obtain  low  p-type  doping  levels  in  the  7i-region  in  HgCdTe 
material  grown  by  molecular  beam  epitaxy  (MBE).  Such  a  low  level  of  p-type  doping  has  never 
been  demonstrated  in  a  controlled  manner  by  MBE  for  HgCdTe  (9-11).  The  goal  of  obtaining  a 
reduced  p-type  carrier  concentration  will  be  to  reduce  the  Auger  generation-recombination 
carrier  dynamics  and  to  provide  the  enhanced  carrier  diffusion  required  for  high  quantum 
efficiency  and  low  dark  currents,  resulting  in  enhanced  D*. 

The  typical  device  structure  is  shown  in  figure  6.  The  intrinsic  carrier  concentration  in 
Hgo.78Cdo.22Te  is  increases  as  the  temperature  is  increased  (see  figure  5).  It  can  reach  values  as 
high  as  1016  cm  3  near  ambient  temperatures,  i.e.  values  higher  than  the  actual  extrinsic  doping 
level.  Therefore,  since  Auger  processes  are  directly  proportional  to  the  number  of  electrons  and 
holes  available,  it  will  result  in  high  dark  currents  in  HgCdTe  IR  devices  operated  at  high 
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temperatures  where  the  absorber  layer  is  intrinsic  at  equilibrium.  The  P+/ti  heterojunction  is 
referred  to  as  the  exclusion  junction  (see  figure  6).  The  j t/  N+  heterojunction  is  referred  to  as  the 
extraction  junction.  For  example,  consider  a  P+/tt/N+  structure  operated  at  high  temperature  so 
that  the  ‘n’  absorber  region  is  intrinsic  in  thermal  equilibrium.  As  an  example  at  200  K,  the 
intrinsic  carrier  concentration  is  in  the  mid-10  cm  for  Hgo.78Cdo.22Te  material  (see  table  1). 
Applying  an  increasing  reverse  bias  to  the  device,  the  7t/N+  extraction  junction  extracts  the 
minority  electrons  from  the  ‘71’  absorber  region  faster  than  they  can  be  injected  from  the  P'/ti 
interface  region.  One  of  the  conditions  is  that  the  diffusion  length  of  electrons  in  the  P+  region 
should  be  less  than  the  thickness  of  P+  region,  and  the  diffusion  length  of  the  electrons  in  the  ‘ n ’ 
absorber  region  should  be  greater  than  the  thickness  of  the  ‘n’  absorber  region.  This  condition 
will  not  allow  electrons  in  the  P+  region  to  cross  over  to  the  ‘n’  absorber  region.  Since  in  the  P+ 
region  electron  concentration,  diffusion  length  of  electrons  is  very  low  (see  the  table),  electrons 
cannot  be  replenished  to  the  ‘71’  absorber  region.  Therefore,  the  electron  concentration  drops 
significantly,  below  its  thermal  equilibrium  value,  in  the  ‘n’  absorber  region.  Higher  band  gap  in 
the  P+  region  help  this  process,  as  well  as  an  electron  barrier  in  between  P+  region  and  ‘71’ 
absorber  region.  The  hole  concentration  is  also  reduced  below  the  equilibrium  value  in  order  to 
maintain  charge  neutrality.  As  the  reverse  bias  is  increased,  the  electron  concentration  in  the 
absorber  region  will  eventually  go  from  the  intrinsic  level  (mid-10  cm  )  down  to  the  extrinsic 
doping  levels,  as  shown  in  figure  5  and  table  1.  This  decrease  of  electron  concentration  in  the 
absorber  region  by  several  orders  of  magnitude  as  the  reverse  bias  is  increased  results  in 
significant  reduction  of  Auger  processes  and,  hence,  reduction  in  dark  current  at  ambient 
temperatures.  Figure  7  illustrates  calculated  electron  and  hole  concentration  in  the  three  regions 
(see  figure  6)  at  300  K. 

Auger  suppression  in  HOT  HgCdTe  devices  translates  into  a  unique  negative  differential 
resistance  (NDR)  in  the  reverse  bias  current-voltage  (TV)  characteristics,  which  has  been 
experimentally  observed  in  devices  grown  by  MOVPE  (6)  and  MBE  using  silver  as  an  acceptor 
dopant  (12).  More  recently,  Wijewamasuriya  et  al  (2008)  reported  Auger  suppression  in  LWIR 
MBE-grown  devices  using  arsenic  as  an  acceptor  dopant  (13). 
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Figure  5.  Temperature  dependence  of  intrinsic 

(solid  line),  electron  (dash  line)  and  hole 
(dotted  line)  concentration  versus  temperature 
of  n-type  Hgo.7sCdo.22Te  material. 
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Figure  6.  HgCdTe  HOT  infrared  photodiode  structure;  “t”  is 
the  P+ layer  thickness,  “d”  is  the  7t-layer  thickness, 
and  “L”  is  the  diffusion  length  of  electrons  in  the 
respective  P+  and  Jt-layers. 
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Table  1.  Calculated  intrinsic,  hole,  and  electron  concentration  in  the  individual  layers  of  P+/7i/N+ 
structure  at  78  K  and  at  200  K. 


Layer 

P 

n 

N+ 

P 

n 

N+ 

X 

0.3 

0.225 

0.3 

0.3 

0.225 

0.3 

at  78K 

at  200K 

nj  (cm3) 

1.50E+09 

7.00E+12 

1.50E+09 

1.00E+14 

3.50E+15 

1.00E+14 

Po  (cm3) 

1.00E+17 

1.00E+15 

2.25E+01 

1.00E+17 

4.00E+15 

4.00E+11 

no  (cm3) 

2.25E+01 

5.00E+10 

1.00E+17 

4.00E+11 

3.00E+15 

1.00E+17 

Figure  7.  Calculated  electron  and  hole  concentration  in  the 

three  layers  (see  figure  6)  of  a  LWIR  HOT  structure 
at  1  V  reverse  bias.  Also  indicated  is  the  doping  levels. 


3.  Carrier  Recombination  Mechanisms  in  HgCdTe 


As  described  earlier,  one  way  to  reduce  dark  current  is  to  increase  the  minority  carrier  lifetime. 
Let  us  look  at  the  carrier  lifetime  mechanisms  in  HgCdTe.  Minority  carrier  lifetime  in  n-type 
HgCdTe  is  determined  by  several  different  recombination  mechanisms,  discussed  extensively  in 
the  literature  (14-27).  In  particular,  Auger,  radiative,  and  Shockley-Read-Hall  (SRH)  processes 
are  important  recombination  mechanisms  in  HgCdTe.  The  first  two — Auger  and  radiative 
mechanisms — are  fundamental  to  the  material,  while  the  last  is  initiated  by  defects  that  may  be 
eliminated  through  high  quality  material  growth  and  device  processing  advances. 
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The  Auger- 1  process  (24)  is  the  dominant  recombination  mechanism  in  n-type  material.  In  this 
mechanism,  an  excited  electron  in  the  conduction  band  recombines  with  a  hole  in  the  valence 
band,  and  the  loss  in  energy  is  transferred  to  a  second  conduction-band  electron  by  the  electron- 
electron  coulomb  interaction.  Overall,  this  involves  two  electrons  and  a  heavy  hole.  Under  low 
injection  conditions,  lifetime  for  Auger- 1  limited  recombination  has  been  described  in  terms  of 

(14), 


K  +  Po)(no  +  4  Po) 

where,  n,  is  the  intrinsic  carrier  concentration,  no,  po  are  equilibrium  electron,  hole  concentration, 
and  zx,  is  the  Auger  lifetime  for  intrinsic  HgCdTe  and  5,  is  the  hole-hole  collision  term,  which  we 
have  neglected  for  HgCdTe.  zx,  is  given  by  (14-16,  19) 


3.8xlQ~18g2(l  +  ^)1/2(l  +  /7)  f  (I  +  27)  £g" 

(melm0)\FxF2\\kTIEgf12  ^  A  (1  +  7)  kT ^ 


giving  lAiin  seconds.  In  this  last  expression,  k  is  Boltzmann’s  constant  and  T  is  the  absolute 
temperature,  Eg  is  the  gap  energy  (28),  and  s  is  the  dielectric  constant,  given  in  (15)  to  be  20  - 
9.4x.  Also,  mo  is  the  electron  rest  mass,  m*e  and  m  /,  are  electron  and  hole-effective  masses,  and 
7  =  m *//«*/,.  The  value  for  IF//7?!  typically  varies  between  0.1  and  0.3  (15-23).  The  best  way  is 
to  get  IF/F?I  empirically  is  to  study  the  concentration  dependence  of  lifetimes  at  78  K  since  there 
exits  ambiguity  in  the  exact  value  of  I F / F? I  in  the  literature.  Figure  8  shows  the  lifetime  behavior 
versus  electron  concentration  of  n-type  doped  HgCdTe  at  78  K  for  X=0.22.  The  lifetime  data 
follows  extremely  well  for  doping  levels  above  10  cm  for  IF;F?I  =  0.18. 

In  p-type  material  in  addition  to  the  Auger- 1  process,  the  Auger-7  process  (24,  28)  is  also 
important.  In  this  process,  an  excited  electron  in  the  conduction  band  recombines  with  a  hole  in 
the  valance  band  and  the  loss  in  energy  is  transferred  to  an  electron  in  the  light  hole  band.  Under 
low  injection  conditions,  Auger-7  lifetime  in  p-type  material  is  given  by  (19,  22), 

T  _  2n?r*M 
K  +  Po )  Po  ’ 

where,  the  ratio  y  =  zA7/  zA1  has  been  calculated  in  (28). 
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Figure  8.  The  dependence  of  minority  carrier  lifetime  on  electron  density  at  78  K 

on  indium-doped  MBE-grown  HgCdTe  layers.  Lines  show  varies  contributions; 
dashed-line  -  SR  recombination  and  red  solid  line  for  p-type  HgCdTe  which  shows 
much  higher  expected  electron  lifetime. 


At  this  time,  there  is  no  clear  agreement  for  y  in  HgCdTe  which  is  related  to  the  magnitude  of 
Auger-7  processes.  Consequently,  we  set  y=  60  in  this  study  which  is  the  most  common  value 
used. 


Next,  radiative  recombination  lifetime  has  been  described  in  terms  of  (15,  16,  19) 

1 

B(n0  +  p0)’ 


with  B  being  given  by 


B  =  5.8x10 


-13£U2 


m« 


V  A 

1  _|_  He.  _|_  He. 


\me  +mh) 


v 


m , 


m 


300 


i 

(e*  +  3kTEg  +3.15k2T2), 


hj 
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giving  tr  i n  seconds.  The  recombination  of  electron-hole  pairs  occurs  between  conduction  and 
valence  bands,  and  the  excess  energy  is  released  in  the  form  of  photons. 


Finally,  Shockley-Read  lifetime,  assuming  a  single  recombination  level  Et  below  the  conduction 
band  edge,  and  low  injection  levels,  has  been  described  in  terms  of  (14,  19,  27), 

vK+^i)  ,  Tn0(Po  +  Pl) 

tsr= - + - ; - > 

n0  +  pQ  n0  +  pQ 


where,  zno  and  zpo  denote  shortest-time  constant  values  for  electron  and  hole  capture, 
respectively,  rij  =  no  exp  [(Et  -  EF)/kT],  pi  =po  exp[(EF-Et)/kT],  and  EF  denotes  the  Fermi 
energy. 

As  can  be  seen  from  figure  8,  for  low  level  of  doping  concentrations  (less  than  10  cm  )  in  n- 
type  material,  SRH  lifetime  is  necessary  in  addition  to  band-to-band  recombination  to  explain  the 
data.  Data  explain  very  well  with  zno  =  zPo  =  5  psec  for  LW  n-type  materials  doped  with  indium. 
This  agrees  very  well  with  data  obtained  in  LPE  materials,  as  well  (2). 

When  the  lifetime  is  not  dominated  by  a  single  process,  it  is  necessary  to  calculate  the  lifetime 
accordingly.  Thus,  if  all  three  recombination  mechanisms  are  important,  the  overall  effective 
lifetime  x  is  then  calculated  from, 

1111 

—  — - 1 - 1 - ^ 

T  T  A  Tr  Tsr 

where  ta,  Tr,  and  tsr  are  the  Auger,  radiative  and  SRH  lifetimes,  respectively. 


4.  Current  Mechanisms  in  a  Photo  Voltaic  Junction 


The  choice  of  PV  detector  architecture  determines  the  type  of  dark  currents  that  impact  detector 
performance,  and  ultimately,  determine  the  signal-to-noise  ratio.  Currents  from  the  n-side  active 
layer  should  only  limit  high-quality  P+/n  devices  operating  in  the  diffusion-limited  region.  In  a 
proper  design  of  a  detector,  the  active  layer  thickness  d  is  less  than  the  diffusion  length  of  the 
minority  carriers.  Dark  diffusion  current  of  a  PV  detector  can  be  expressed  as  follows  (29): 

lDiff  =  W 


exp 


qv 

nkT 
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where  q  is  the  electron  charge,  k  is  Boltzmann’s  constant,  T  is  absolute  temperature,  n  is  the 
ideality  factor  equals  to  1,  Vis  the  bias  voltage  across  the  diode,  and  JoDiff  is  the  saturation 
diffusion  current  density. 


Diffusion  current  contributions  to  reverse  saturation  current  JoDijf  consist  of  components  from  p 
and  n-sides  of  the  junction,  expressed  here  in  current  density  superposition  form  as  (30-33) 


J, 


Odiff 


=  Vh 


L. 


Ndrp 


+ 


L. 


where  J  odiff  denotes  diffusion  current  density,  Lpn  and  tp>„  denote  minority  diffusion  lengths  and 
lifetimes  on  respective  n  and  p  sides  of  the  junction  ,  and  Nd,a  denote  donor  and  acceptor 
concentration  levels  on  the  n  and  p  sides  of  the  junction.  For  the  case  of  p-on-n  double 
heterostructures  where  the  p-side  has  a  wider  bandgap  and  high  doping  N+a ,  the  second  term  in 
Jo  diff  is  negligible  compared  to  the  first  and  J<knff  reduces  to 

cin!  LP 

0diff  Ndr 

d  p 

When  the  n-type  active  (absorber)  layer  thickness  d  is  smaller  than  the  diffusion  length  Lp  , 
dimension  d  effectively  replaces  Lp  to  give 

qnfd 
J  0diff  =  -jj— 

dl  P 

Second  important  current  mechanism  is  the  generation-recombination  (g-r)  current  density  which 
arises  due  to  imperfections  in  the  space  charge  region. 

Early  treatments  developed  the  seminal  formulation  for  g-r  phenomena  in  semiconductor 
junctions  (34).  Phenomenologically,  electron-hole  pairs  generated  thermally  in  the  depletion 
region  give  rise  to  a  current  that  is  recombination-behaving  under  forward  bias  and  generation¬ 
behaving  under  reverse  bias.  Derivations  for  this  current  historically  take  the  form; 


J 


g-r 


qnyV 


J 


TnOrpO  q 


K-v), 


2  kT 


fib) 


where, 
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/M=J 


du 


u  +  2  bu  + 1 


b  =  e 


-qV, 


'2kT 


cosh 


E,~Ei  1 

— - -  +  —In 

kT  2 


1  pO 

V  Tn0  J 


Trap  levels  due  to  defects  are  most  effective  when  they  exist  close  to  the  intrinsic  Fermi  level, 
i.e.,  when  E,  =  E,.  For  bias  conditions  where  V  ~0,  along  with  assumptions  that  E,  =  E,  and  zpo  = 
zno  =  zm  it  follows  that  b  »1  and  f(b)  ~  /.  Under  these  provisos,  Jg.r  reduces  to 


7  = 


qnyV  2  kT 


sinhl 


/  \  JLJ  1  1  1  I 

q{vb,-v)  \ 


K 


2  kT 


Generation-recombination  current  Jg.r  current  is  seen  to  be  proportional  to  n,.  As  such, 
generation-recombination  J,,.r  current  will  tend  to  dominate  over  diffusion  /^current  at  low 
temperatures. 

From  the  previous  equation,  g-r  current  arises  due  to  imperfections  in  the  space  charge  region 
and  takes  the  form  of 


I  g-r  I  og-r 


exp 


qV 

nkT 


-1 


where  Iog.r  is  the  g-r  saturation  current  (which  is  a  fitting  parameter)  and  n  is  the  ideality  factor 
equal  to  2. 


Tunneling  currents  can  dominate  in  narrow-band  gap  HgCdTe  devices  at  low  temperatures, 
where  thermally  generated  currents  are  low,  and  this  is  the  dominant  current  mechanism  under 
reverse  biases  under  such  conditions.  Tunneling  currents  come  in  two  forms,  namely,  (1)  band- 
to-band  (BTB)  and  (2)  trap-assisted  tunneling  (TAT),  where  forbidden  gap  energy  states  promote 
and  participate  in  a  two-step  process.  BTB  tunneling  currents  result  in  electrons  tunneling 
directly  from  the  valance  band  to  the  conduction  band  and  are  given  by  (55) 


Aq  EV  2m  ± 
BTB  4x2H2  V  E„ 


exp 


4x/2 


m  e: 


3qEh 


1-exp 


'  qV 
KkT 


where  A  is  the  junction  area,  E  is  the  electric  field,  ti  is  the  reduced  Planck  constant,  m*  is  the 
tunneling  effective  mass,  and  Eg  is  the  band  gap.  These  currents  are  due  to  direct  tunnelling  of 
carriers  under  relatively  high  reverse  biases.  These  currents  are  fairly  independent  of 
temperature,  and  calculated  BTB  currents  at  78  K  for  a  band  gap  of  1 15  meV  is  insignificant. 

TAT  current  (Itat)  arises  from  electrons  transitioning  from  the  valence  band  to  a  trap  level 
within  the  forbidden  gap  and  then  tunneling  to  the  conduction  band.  The  current  density  is  given 
by  (36) 


13 


Itat  — 


Art2q2m  EM  2NtW  2 m 


n\E-E, ) 


exp 


4m  12  Elg5F(a ) 

2qEh 


1-exp 


qV 

kT 


where  M  is  the  transition  matrix  element,  N,  is  the  trap  density  in  cm  ,  Et  is  the  trap  energy,  and 
F(a)  is  function  of  E/E,,. 

Finally,  currents  due  to  shunting  (Ish)  across  the  diode  are  given  by 


where  Rsh  is  the  shunt  resistance.  Usually,  surface  leakage  currents  and  shunt  currents  are  due  to 
the  intersection  of  dislocation,  defects,  etc.,  which  are  possible  sources  of  this  category  of  excess 
currents. 


Combining  the  previously  mentioned  considerations  relative  to  diffusion,  g-r,  BTB,  TAT,  and 
shunt  currents,  detector  overall  dark  current  is  given  by 

I  Total  ~  I D  iff  I  g-r  ^  BTB  ^  TAT  ^  SH  ■ 


Figures  9  and  10  compare  the  measured  versus  modeled  I- Vs  from  the  previous  equations  for 
two  P+/n  HgCdTe  PV  devices.  The  model  calculates  ideal  diffusion  currents  from  the  active 
region  from  the  measured  material  parameters,  such  as  doping  Nd,  active  layer  thickness  d,  and 
Ac  at  78  K.  The  measured  Ac  at  78  K  is  10.0  pm.  Diffusion  currents  from  the  n-type  active  region 
fit  very  well  down  to  -0.05  V,  indicating  that  the  dominant  current  mechanism  is  diffusion  close 
to  zero  bias  and  a  diffusion  current  of  3. 18  x  10  10A  was  measured.  At  larger  biases  (greater 
than  -0.4  V),  the  measured  currents  were  significantly  higher  than  the  diffusion  currents,  and  our 
model  calculation  shows  that  these  are  due  to  TAT  (dotted  line).  The  solid  line  in  figure  9  is  the 
contribution  of  diffusion,  TAT,  and  shunt  currents.  The  BTB  and  g-r  current  contributions  were 
negligible.  A  very  good  agreement  exists  between  the  measured  and  model  calculations,  as  can 
be  seen  from  figure  9.  As  evidenced  by  the  Rd  versus  V  curve,  there  is  a  slight  shunt  current 
contribution  between  the  -300  and  -50  mV  region.  Shunt  resistance  is 
5xl09  Q.  Close  to  zero  bias,  the  obtained  shunt  resistance  is  much  higher  than  the  diode 
impedance  and,  hence,  has  no  effect  on  zero  bias  impedance  (RqA)  on  “good  devices.”  On  the 
other  hand,  figure  10  compares  the  measured  versus  calculated  I-Vs  for  another  diode.  In 
contrast  to  the  previous  diode,  a  flat  reverse  saturation  current  is  not  seen  close  to  zero  bias.  The 
reverse  current  linearly  increases  as  the  bias  voltage  increases,  and  ultimately,  at  larger  biases, 
TAT  currents  dominate.  In  order  to  get  a  good  fit  to  the  measured  I-V  data,  a  shunt  current  plus 
TAT  current  are  needed  in  addition  to  the  diffusion  current.  Currents  are  diffusion  at  zero  bias, 
but  a  slight  reverse  bias  shows  a  deviation  from  this  current  mechanism.  Instead,  measured 
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currents  close  to,  but  not  at,  zero  biases  exhibit  a  shunt  behavior  and  shunt  resistance  of 
7.8  107  Q,  which  have  been  assumed  in  order  to  get  a  good  fit  to  the  measured  data. 


Figure  9.  Dependence  of  measured  (symbols)  and  modeled  (lines)  of  current  versus  bias  voltage  at 
78  K.  The  total  current  (solid  line)  is  modeled  with  a  diffusion  saturation  current  of 
3.18x10  10  Amps  along  with  trap-assisted  tunneling  and  shunt  currents.  A  shunt  resistance 
of  5xl09  Q  is  assumed  in  the  model  calculation. 

The  dashed  line  in  figure  10  with  n  =  1  matches  the  R&  versus  voltage  curve,  further  confirming 
that  the  dominant  current  mechanism  is  diffusion. 

The  lowest  dark  current  is  achieved  when  radiative  recombination  is  dominant.  The  radiative 
recombination  process  has  been  found  to  be  a  non-issue  in  near-equilibrium  situations  and  can  be 
ignored  due  to  photon  recycling  (2,  37).  This  leaves  generation  currents  from  the  SRH  centers  as 
the  major  source  of  leakage  currents  limiting  performance,  even  at  ambient  temperatures.  Hence, 
the  ultimate  limitation  is  the  Auger  recombination  and,  if  suppressed,  this  will  lead  closer  to 
background-limited  directivity  D*  at  ambient  temperatures  on  good-quality  p-type  absorber 
materials.  Figure  11  shows  the  calculated  saturation  dark  current  density  versus  temperature  for 
Xco  =  10  pm  photovoltaic  devices  at  100  K  for  a  doping  density  of  2.0x10  cm  .  The  Cd 
composition  is  0.22.  The  PV  device  designs  under  consideration  are  double-layer  hetero structure 
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P+/n,  N+/p,  and  PVti/N+  under  high  reverse-bias  operation.  As  shown  in  the  figure,  high  dark 
current  dominates  at  operation  near  ambient  temperature  in  P+/n  structures  due  to  high  Auger 
recombination  in  n-type  materials.  This  dark  current  could  be  decreased  in  N+/p  devices  due  to 
higher  minority  Auger  carrier  lifetimes  (factor  of  y)  than  n-type  absorber  layers.  This  dark 
current  could  be  suppressed  even  further  by  the  phenomena  of  exclusion  and  extraction  in 
P+/ti/N+  devices,  as  indicated  in  the  figure.  The  phenomena  of  minority  carrier  exclusion  and 
extraction  under  high  biases  have  been  described  previously  {1 ,  6-8).  In  the  P+/tt/N+  device 
architecture,  the  jr-region  is  the  absorbing  layer.  At  ambient  temperatures,  the  minority  and 
majority  carrier  concentrations  in  the  p-region  are  approximately  equal  and,  hence,  behave  as 
near  intrinsic  (see  table  1).  This  will  give  high  dark  current,  and,  therefore,  noise,  due  to 
dominating  intrinsic  carriers. 


Figure  10.  Dependence  of  measured  (symbols)  and  modeled  (lines)  of  current  versus  bias  voltage  of  a 
diode  at  78  K.  The  total  current  (solid  line)  is  modeled  with  diffusion  saturation  current  of 
1.08  x  10  10  Amps  along  with  trap  assisted  tunneling  and  shunts  currents.  A  shunt 
resistance  of  7.8  x  107  O  is  assumed. 

Trick  is  to  reduce  the  minority  carrier  concentration — as  an  example  at  200  K,  from  3xl016  cm"3 

1 S  _'2 

to  concentration  levels  of  3x10  cm  at  78  K. 

The  concept  is  to  reduce  both  minority  and  majority  carrier  concentrations  from  their  equilibrium 
value  to  the  levels  of  extrinsic  conditions  so  that  extrinsic  behavior  will  be  observed  at  ambient 
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temperatures  (6-8).  These  reduced  carriers  will  give  suppressed  Auger  recombination  lifetimes. 
Figure  1 1  shows  the  dark  current  density  versus  temperature  on  such  a  P+/ti/N+  device 
architecture  under  exclusion  and  extraction.  As  can  be  seen,  this  device  architecture  gives  the 
lowest  dark  current  densities  as  ambient  temperatures  are  approached,  which  gives  the  potential 
to  yield  closer  to  BLIP  operation.  Further,  if  very  low  levels  of  p-type  acceptor  concentrations 
(Na  ~  1014  cm-3)  can  be  achieved,  then  near-BLIP  operation  at  ambient  temperature  should  be 
possible.  The  photocurrent  due  to  300  K  blackbody  radiation  is  4xl0-5  A  for  a  quantum 
efficiency  of  unity.  As  Elliott  and  co-workers  have  pointed  out,  one  of  the  main  limitations  of 
this  structure  is  the  leakage  of  electrons  from  the  P+  to  the  p-region  close  to  ambient 
temperatures.  A  slightly  modified  structure  has  been  proposed  to  surmount  this  situation  (38).  In 
this  report,  we  discuss  the  results  of  P+/ti/N+  device  structures  grown  by  MBE. 
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Figure  11.  Calculated  dark  current  density  versus  temperature:  (a)  for  P+/n 

device  (dotted  line),  (b)  for  N+/p  (dash  line),  and  (c)  P+/7i/N+  (solid 
line)  under  high  reverse-bias  operation.  The  calculation  assumes 
a  one-dimensional  (1-D)  device  model,  currents  only  from  active 
absorber  layers,  and  active  absorber  doping  concentrations  of 
2xl015  cm  3.  For  comparison  purposes,  the  photocurrent  due  to 
300  K  blackbody  radiation  is  4x  1 0  3  Amps  for  a  quantum  efficiency 
of  unity. 
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5.  Experiment 


5.1  MBE  Growth  of  HgCdTe 

Layers  characterized  here  were  grown  by  MBE  on  (211)B  CdZnTe  substrates.  Figure  12  shows  a 
schematic  cross  section  of  the  P'/ti/N+  planar  detector  architecture.  The  key  feature  of  the 
detector  architecture  is  the  planar  device  geometry,  with  a  narrow-bandgap  active  p-layer 
sandwiched  between  a  wide-band  P+  cap  layer  and  a  wide-band  N+  layer.  Benefits  of  this 
architecture  include  a  buried  narrow-bandgap  active  layer,  resulting  in  a  reduction  in  surface 
generation-recombination  and  surface  tunneling  currents,  and  an  inherent  ease  to  passivate  the 
wider-band  gap  cap  layer  (39,  40).  ^indicates  higher  Cd  composition  as  well  as  higher  p-type 
doping  than  the  base  layer  and  N+  indicates  higher  Cd  composition,  as  well  as  higher  n-type 
doping  than  the  base  layer.  All  three  regions  shown  in  figure  12  are  in-situ  indium-doped  during 
the  MBE  growth,  and  the  doping  levels  are  shown  in  the  table  in  the  figure.  Arsenic-ion 
implantation  was  performed  in  selected  regions  followed  by  a  deep  diffusion  plus  activation 
anneal  to  fabricate  P+/tt/N+  detectors  as  a  planar  configuration  into  process  evaluation  chips 
(PEC)  containing  devices  with  variable  areas. 


Figure  12.  (a)  Cross  section  of  the  P+/7t/N+ planar  device  architecture  grown  by  molecular  beam 
epitaxy.  Implantation  plus  a  diffusion  process  were  used  to  obtain  p-absorber  layer. 

(b)  Table  of  material  characteristics  of  the  different  regions  of  the  device  structure  before 
and  after  implantation/diffusion  anneal. 
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We  mainly  concentrated  on  250-pm-diameter  devices  in  this  report.  These  fabricated  devices 
were  then  characterized  electrically  and  optically  down  to  78  K  from  300  K.  Backside 
illumination  was  used  for  all  optical  measurements.  The  spectral  response  was  measured  with  a 
Fourier  transform  infrared  (FTIR)  spectrometer,  and  the  50%  spectral  response  point  was  taken 
as  the  cut-off  wavelength  of  the  device. 

5.2  Low  Arsenic  Doping  in  HgCdTe  -  n  Absorber  Region 

In  this  work,  arsenic  was  used  as  a  source  of  extrinsic  p-type  carriers  with  the  proper  arsenic 
doping  level  obtained  by  a  novel  diffusion  technique  based  on  diffusion  by  thermal  annealing  to 
obtain  PVti  region  in  the  planar  devices  architecture  in  MBE-grown  n-type  HgCdTe  layers. 

In  early  work  (72),  Auger  suppression  in  MBE-grown  HgCdTe  heterostructure  devices  using 
silver  (Ag)  as  a  p-type  doping  was  demonstrated.  However,  the  use  of  group  V  elements  as  p- 
type  dopants  in  MBE  has  been  studied  extensively  due  to  the  larger  atomic  size  of  the  group  V 
elements  and,  hence,  their  lower  diffusivity  in  HgCdTe.  This  fact  allows  for  stable  and  well- 
controlled  p-n  junctions  with  sharp  transition  interfaces  (41).  The  necessary  control  of  low 
1015  cm  1  p-type  doping  for  the  P'/ti:/N+  structure  has  never  been  achieved  in  more  than  20  years 
of  research  on  HgCdTe  development  by  MBE.  Present  technology  limits  p-type  doping  with  As 
values  down  to  ~1016  cm" 3  and  even  this  amount  of  doping  is  usually  not  obtained  reproducibly. 
Our  novel  diffusion  process  allowed  arsenic  concentrations  well  below  the  present  technology. 

Figure  13  shows  a  Gaussian-like  distribution  of  an  arsenic-implanted  arsenic  (open  circles) 
reaching  a  depth  of  0.6  pm  from  the  surface.  Deep  diffusion  of  arsenic  is  achieved  after 
annealing  under  a  very  low  partial  pressure  of  Hg,  as  initially  demonstrated  by  Wijewarnasuriya 
et  al.  (13).  As  an  example,  figure  13  shows  the  arsenic  profile  by  secondary  ion  mass 
spectroscopy  (SIMS)  after  16  h  of  annealing  at  300  °C  under  appropriate  Hg  partial  pressure 
(solid  circles).  After  deep  diffusion,  the  arsenic  depth  reaches  5  pm.  After  the  diffusion  anneal, 
an  activation  anneal  was  conducted  at  400  °C  for  5  min  under  Hg  saturated  conditions,  followed 
by  a  Hg  vacancy  annihilation  anneal  at  250  °C  for  18-20  h  under  Hg  saturated  conditions.  The 
alloy  composition  was  not  measured  by  SIMS,  and  the  author  is,  therefore,  not  able  to  comment 
on  Cd  interdiffusion.  An  important  condition  of  P+ region  is  that  minority  carrier  concentration 
should  be  low;  consequently,  under  reverse  biases,  electrons  should  not  replenish  from  the  P+ 
region  to  the  p-absorbing  layer. 
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Figure  13.  Typical  example  of  arsenic  redistribution  from  an 
arsenic-ion  implanted  source.  Open  symbols  are 
as-implanted  arsenic  measured  by  SIMS  and  solid 
symbols  show  arsenic  concentration  after  annealing 
at  300  °C  under  Hg  diffusion  conditions  for  16  hours. 

Lines  are  Gaussian  fits  to  the  data. 

The  main  parameters  of  variation  are  the  anneal  duration,  Hg  source  temperature  and 
corresponding  Hg  partial  pressure  P(Hg).  All  samples  were  held  at  300  °C  during  the  deep 
diffusion  anneal.  For  each  condition,  the  samples  were  sealed  in  a  quartz  ampoule  and  placed 
into  a  two-zone  annealing  apparatus,  one  zone  for  the  samples  and  the  other  zone  for  the  Hg 
source.  The  partial  pressure  of  Hg  was  controlled  by  the  Hg  source  temperature,  and  the 
corresponding  pressure  values  were  determined  from  a  look-up  table  (42).  Activation  and 
vacancy  annihilation  anneals  were  carried  out  at  400  °C  for  5  min  and  at  250  °C  for  20  h, 
respectively.  A  5  °C  difference  between  the  sample  and  Hg  source  was  maintained  to  prevent 
condensation  of  Hg  droplets  on  the  sample  surface  during  quenching.  The  As-diffusion  profiles 
in  the  HgCdTe  layers  were  measured  by  SIMS  at  Evans  Analytical  Group.  As  diffusion 
coefficients  D(As)  were  then  extracted  from  the  profiles  using  a  dopant  distribution  model,  as 
described  below  (43): 


(  (d-R)2^ 

G  I -  EAjJ 

■dnDt 

{  4  Dt  J 

where,  d  is  the  depth  from  the  sample  surface,  t  is  the  annealing  duration,  Q  is  the  implant  dose, 
and  R  is  the  fitting  parameters  for  peak  depth. 


20 


The  SIMS  profiles  for  LWIR  material  annealed  with  varying  Hg  partial  pressure  are  shown  in 
figure  14. 


Figure  14.  SIMS  data  of  LWIR  As  diffusion  profile  for  (a)  16  hour  anneal  at  P(Hg)  0.07  atm  and  (b)  20.5  hour 
anneal  at  P(Hg)  =  0.32  atm. 

Comparing  the  junction  depth  of  the  P+/7T  interface  and  the  absorber  region  widths,  deeper  As 
diffusion  is  observed  in  the  sample  annealed  under  lower  Hg  pressure  (0.07  atm).  The  data  in 
figure  14a  exhibits  a  junction  depth  of  2.5  pm,  an  absorber  region  between  the  depths  of  1  pm 
and  5  pm,  and  an  As  doping  concentration  in  the  range  of  1015  -  10 16  cm"  \  The  data  in  figure 
14b  exhibits  a  junction  depth  of  0.6  pm,  an  absorber  region  between  0.8  pm  and  1.6  pm,  and  an 
As  doping  level  in  the  range  of  1016  -  1017  cm  3.  The  As  concentration  in  the  absorber  region  is 
an  order  of  magnitude  lower  for  the  sample  annealed  at  a  lower  Hg  pressure  due  to  the  increased 
diffusion  of  As  atoms  through  a  higher  concentration  of  Hg  vacancies.  As  diffusion  coefficients 
of  D(As)  =  6xl0-14  cm2S  1  and  D(As)  =  2xl0-15  cm2S  1  are  extracted  for  samples  in  figures  14a 
and  14b,  respectively, 

All  samples  exhibited  consistent  P+/7t  step-junction  profiles  with  an  arsenic  concentration 
contrast  of  approximately  two  orders  of  magnitude  in  the  P+  and  n  regions.  We  observe  two 
primary  components  in  the  arsenic  diffusion  profiles:  (1)  a  Gaussian  component  and  (2)  a  linear 
component,  as  indicated  in  figure  15.  These  characteristics  contribute  to  the  P+/7i  profile  in  which 
we  are  interested,  As  reported  in  past  studies  (44),  the  movement  of  As  in  HgCdTe  during 
annealing  can  be  described  by  fast  and  slow  diffusion  components.  We  believe  the  Gaussian 
distribution  profile  is  attributed  to  the  slower  diffusion  of  arsenic  atoms  traveling  primarily  via 
Hg  vacancies.  We  also  believe  that  the  linear  tail  of  the  profile  may  be  a  result  of  arsenic 
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channeling  and  possibly  diffusion  by  means  of  other  local  defects  associated  with  the  ion 
implantation  process. 


Figure  15.  SIMS  data  of  LWIR  As  diffusion  profiles  of  as-implanted,  after  two  hours,  and  after  four  hours  of 
annealing  at  300  °C  under  Hg-deficient  condition.  Obtained  D(As)  =  9x10  13  crn2S 

Our  diffusion  coefficients  follow  the  general  trend  of  the  inverse  relationship  between  the  arsenic 
diffusion  coefficient  and  the  Hg  partial  pressure  as  described  in  past  diffusion  experiments  (45— 
47).  Figure  16  compares  our  extracted  diffusion  coefficients  against  values  obtained  in  past 
studies. 
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Figure  16.  Published  arsenic  diffusion  coefficient  in  HgCdTe. 

It  is  important  to  note  that  our  data  points  correspond  to  MBE-grown  LWIR  samples  annealed  at 
a  low  substrate  temperature  of  300  °C  and  longer  annealing  times,  while  most  of  the  past  studies 
focused  on  LWIR  samples  grown  by  MOCVD  and  LPE  methods  and  annealed  for  short 
durations  at  higher  substrate  temperatures.  We  notice  there  is  a  difference  in  the  diffusion 
coefficients  we  obtained  for  the  LWIR  samples.  Comparing  our  LWIR  diffusion  coefficients  for 
samples  annealed  at  300  °C  to  those  projected  by  the  diffusion  coefficient  model  proposed  for 
Hgo.8Cdo.2Te  by  Shaw  et  al.  (48)  for  the  same  temperature,  we  observe  that  our  data  points  lie 
close  to  the  model.  Slight  differences  in  the  D(As)  values  may  correspond  primarily  to  non¬ 
idealities  introduced  by  the  actual  experiment.  Lurther  studies  need  to  be  done  to  compare  the 
diffusion  coefficients  of  LWIR  samples  annealed  at  lower  substrate  temperatures  to  existing 
diffusion  coefficient  models. 

A  critical  aspect  of  our  novel  diffusion  process  is  the  ability  to  maintain  a  good  HgCdTe  surface 
even  after  successful  deep  diffusion  of  As  using  a  low  Hg  partial  pressure.  Several  experiments 
were  also  carried  out  under  different  Hg  pressures  to  optimize  the  surface  smoothness  of  the 
HgCdTe  surface  before  and  after  annealing.  Ligure  17  shows  the  surface  under  optimized  arsenic 
deep  diffusion  conditions  before  and  after  annealing.  The  HgCdTe  surface  will  become  rough  if 
the  arsenic  diffusion  anneals  conditions  are  not  optimized,  as  shown  in  figure  18. 
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In  order  to  take  advantage  of  the  higher  minority  carrier  lifetime  in  p-type  material,  the  absorber 
layer  has  to  be  lightly  doped.  This  was  achieved  using  our  unique  deep  diffusion  process  at 
300  °C,  as  indicated  in  figures  13  to  15.  The  doping  concentration  gradient  in  region  II  also  helps 
to  further  diffuse  electrons  (in  this  case  minority  carriers  in  region  II)  to  the  p/N+ junction.  The 
thickness  of  the  absorber  layer  chosen  in  this  case  is  about  5  pm,  which  is  less  than  the  diffusion 
length  of  minority  electrons  in  region  II.  This  is  true  even  at  ambient  temperature,  since  the 
mobility  of  minority  carriers  (in  these  case  electrons)  is  on  the  order  of  high  10  cm7V  s  and 
lifetime  values  are  on  the  order  of  microseconds. 


Figure  17.  (a)  Surface  morphology  under  Normaski  microscopy  before  and  after  diffusion  anneals  and 
(b)  showing  implanted  regions  under  Normaski  after  complete  anneals. 
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After  Anneal 


Rough  Surface  at  the 
bottom  of  trenches 


Figure  18.  Surface  morphology  of  the  trenches  under  Normaski  microscopy  before  and  after  diffusion  anneals. 


6.  Current- Voltage  Characteristics 


I — V  measurements  were  made  from  biases  between  -1000  mV  and  0  mV.  These  measurements 
were  made  using  a  HP  parameter  analyzer.  Figure  19  shows  such  I-V  measurements  on  three 
250  pm-diamctcr  devices  at  300  K.  When  the  reverse  bias  increases,  a  sudden  drop  in  the  reverse 
current  is  observed.  This  sudden  drop  is  attributed  to  Auger  suppression  and  is  due  to  the 
removal  of  the  electrons  in  region  II.  As  the  electrons  are  removed  from  region  II,  hole 
concentrations  are  also  reduced  below  the  equilibrium  value  in  order  to  maintain  the  charge 
balance.  Hence,  one  of  the  critical  interfaces  of  this  type  of  architecture  is  the  P+/p  excluding 
interface.  Leakage  currents  associated  with  electron  replenishing  at  this  interface  are  one  of  the 
major  sources  of  high  currents  at  ambient  temperatures.  A  slightly  modified  architecture  has 
been  proposed  to  surmount  this  situation  (38).  At  the  transition  voltage  point,  the  devices  exhibit 
infinite  impedance.  From  the  two  measured  devices,  a  relatively  small  spread  in  currents  can  be 
observed,  indicating  good  uniformity  in  doping  and  material  characteristics.  These  I-V  curves 
show  leakage  currents  with  no  soft  breakdown  characteristics  in  the  high  reverse-bias  region, 
further  indicating  good-quality  MBE  material  and  fabrication  process.  An  approximately  50% 
reduction  in  reverse  current  at  high  biases  is  shown  in  the  figure.  As  indicated  in  the 
introduction,  the  main  limitation  of  this  structure  is  the  leakage  of  electrons  from  the  P+  to  the  p- 
region.  An  approximately  50%  reduction  in  Auger  suppression  is  seen,  using  Ag  as  a  p-type 
dopant  in  a  similar  structure  grown  by  MBE  (12),  and  our  present  work  also  indicates  about  a 
50%  reduction  at  300  K. 

To  assess  the  mechanisms  controlling  the  device,  I-V  characteristics  were  measured  on  device  as 
a  function  of  temperature  down  to  130  K.  As  seen  in  figure  20,  this  device  exhibited  negative 
resistance  characteristics  from  300  K  down  to  130  K.  Figure  19  shows  the  ratio  of  Imax  to  Imin  as  a 
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function  of  temperature.  At  ambient  temperatures,  the  ratio  tends  to  saturate  at  2.  A  similar  trend 
was  observed  in  MOVPE  grown  structures  by  Elliott  and  co-workers  (6-8). 


Figure  19.  Measured  current  voltage  characteristics  at  300  K  showing 
IMax,  IMin  and  negative  resistance  region  which  resulted  from 
suppression  of  Auger  recombination. 

NOTE:  At  high  negative  biases  (away  from  negative  resistance  region)  measured  current  is 
independent  of  applied  bias  indicating  negligible  tunneling  currents.  Calculated  series  resistance 
from  the  slope  closer  to  zero  bias  is  78  Ohms.  Device  area  is  Aj  =  4.9x10  4  cm2.  One  of  the  two 
measured  device  at  300  K  showing  a  steep  transition  from  IMax  to  IMin  in  the  reverse  bias  (solid 
triangular  symbol).  It  is  not  well  understood  why  some  devices  show  a  sharp  transition  and  some 
devices  show  a  soft  transition  from  IMax  to  IMin- 

Spectral  response  versus  wavelength  was  measured  under  backside-illuminated  conditions  on  the 
device.  Figure  21  shows  the  response  versus  wavelength  on  the  device  at  100  K.  The  measured 
cut-off  wavelength,  Xc,  is  approximately  10  pm  at  100  K.  The  calculated  response  per  photon 
versus  wavelength  is  also  displayed  in  figure  21  and  can  be  seen  to  be  near-classical  behavior. 
The  slight  drop  in  response  at  shorter  wavelengths  is  attributed  to  the  interface  recombination  on 
the  N+  layer.  The  detectors  are  non-antireflection  coated  at  the  back  of  the  CdZnTe  substrate  and, 
hence,  approximately  22%  absorption  is  lost  to  reflection  losses.  The  measured  QE  at  -0.1  V  at 
7  pm  wavelength  is  0.54  at  100  K,  which  indicates  the  excellent  quality  of  the  HgCdTe.  Softness 
of  the  response  close  to  the  cut-off  wavelength  region  is  due  to  the  insufficient  absorption  due  to 
the  short  absorber  layer  thickness. 
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At  100  K,  this  device  showed  contact  issues  due  to  less  than  desired  p-type  doping  at  the  top 
contact  (see  figure  13).  This  is  confirmed  by  the  measured  low  QE  at  zero  bias  at  100  K.  We 
have  also  measured  response  at  250  K  (figure  22)  and  a  near  classical  response  was  observed. 
This  indicates  that  diffusion  or  collection  lengths  are  larger  than  the  absorber  layer  thickness.  At 
250  K,  the  diffusion  length  is  expected  to  be  lower  than  at  100  K,  making  it,  therefore,  surprising 
that  there  is  a  signal  even  at  ambient  temperatures  such  as  250  K.  The  measured  cut-off  is 
approximately  7.4  pm  at  250  K  (figure  23). 
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Figure  20.  Measured  ratio  of  Imax  to  I,n[n  as  a  function  of  temperature. 
The  ratio  reaches  ~2  at  300  K. 
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Figure  21.  Measured  current- voltage  versus  temperature  characteristics  of  the  device  (see  text) 
indicating  negative  differential  resistance  down  to  130  K. 
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Figure  22.  Measured  spectral  response  versus  wavelength  at  100  K  for  a  250 -pm- diameter  device. 

Both  spectral  responses,  per  photon  and  per  Watt,  show  near -classical  behavior. 

NOTE:  A  slightly  reduced  response  at  shorter  wavelengths  is  indicative  of  some  surface 
recombination  at  the  back  interface.  The  obtained  cut-off  (A,c)  is  10.4  pm  at  100  K.  The  cut-off  is 
defined  as  the  50%  spectral  response  point  from  the  peak  wavelength.  The  measured  non- 
antireflective-  coated  device  quantum  efficiency  at  -0.1  V  is  0.57  at  a  wavelength  of  7  pm.  The 
measured  quantum  efficiency  at  zero  bias  showed  contact  issues  due  to  a  less  than  desired  p-type 
doping  concentration  at  the  top  contact  (see  figure  13). 
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Figure  23.  Measured  spectral  response  versus  wavelength  at  250  K  for  a  250-pm-diameter  device. 
Both  spectral  responses,  per  photon  and  per  Watt,  show  near-classical  behavior.  The 
soft  transition  close  to  the  cut-off  is  due  to  insufficient  absorption  due  to  the  5  pm 
thin  p-layer.  The  obtained  cut-off  0.c)  is  7.4  pm  at  250  K. 


7.  Fitting  Methodology 


Steady-state  numerical  simulations  are  performed  using  Sentaurus  Device,  a  commercial 
package  by  Synopsys  (49).  Sentaurus  Device  self-consistently  determines  a  full  coupled  solution 
of  the  Poisson’s  equation  and  the  continuity  equations  for  electrons  and  holes  using  an 
approximate  Newton  method  with  Gaussian  elimination.  Fermi-Dirac  statistics  are  used. 

HgCdTe  material  and  recombination  parameters — considering  Auger,  radiative  and  SRH 
recombination — are  included  in  the  model.  Carrier  degeneracy  and  conduction  band  non- 
parabolicity  are  also  taken  into  account.  Surface  recombination  and  tunneling  mechanisms  are 
not  accounted  for  in  the  simulations.  Simulation  output  calculations  of  electrical  and  optical 
characteristics  of  HgCdTe  p-n  infrared  photodiodes  are  comparable  with  published  experimental 
data  from  different  sources  for  long-wave  infrared  (LWIR)  (50).  The  Cd  alloy  profile,  including 
compositional  grading  at  the  two  junctions,  is  taken  from  experimentally  measured  spectroscopic 
ellipsometry  (SE)  data  and  included  in  the  simulation  using  a  linear  interpolation  of  the  HgCdTe 
material  parameters. 
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The  calculated  energy  band  diagram  for  this  device  at  1  V  reverse  bias  and  T=300  K  is  shown  in 
figure  24a.  Partial  depletion  of  the  absorber  layer  is  achieved  by  increasing  the  reverse  bias 
(depletion  width  is  ~2-2.5  pm).  The  calculated  equilibrium  and  non-equilibrium  (IV  reverse 
bias)  carrier  concentration  is  shown  in  figure  24b.  Under  equilibrium,  the  LWIR  absorber  layer  is 
essentially  intrinsic.  As  the  reverse  bias  is  increased,  the  electron  and  hole  concentration 
decreases  by  several  orders  of  magnitude  due  to  the  effect  of  the  exclusion/extraction  layers.  The 
hole  concentration  decreases  down  to  the  extrinsic  p-type  doping  level  for  a  depth  <4  pm  (non- 
depleted  region).  This  initial  simulation  confirms  that  the  NDR  observed  experimentally  is  likely 
due  to  strong  Auger  suppression  in  the  device. 
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Figure  24.  (a)  Planar  device  structure,  (b)  Cd  alloy  composition  after  annealing  obtained  by  SIMS,  and  (c)  doping 
profile  in  the  HOT  structure. 

Our  device  model  is  used  to  fit  the  temperature-dependent  current-voltage  data  of  HOT  #42 
shown  in  figure  19.  We  perform  a  ID  simulation  of  a  cross-section  of  the  device.  In  order  to  fit 
the  I-V  data  for  each  temperature  (120  K,  130  K,  150  K,  200  K,  270  K,  and  300  K),  we  fit  the 
maximum  leakage  current  experimental  value  before  Auger  suppression  (Jmax)  and  the  minimum 
leakage  current  experimental  value  after  Auger  suppression  (Jmin)  as  illustrated  in  figure  19.  In 
our  model,  Jmm  is  taken  as  the  calculated  current  density  at  500  mV  reverse  bias,  which  is  found 
to  be  high  enough  for  dark  current  saturation  (neglecting  the  series  resistance).  Jmax  is  taken  as 
the  maximum  calculated  current  density  between  zero  bias  and  500  mV  reverse  bias.  The 
temperature-dependent  experimental  values  for  Jmin  and  Jmax  that  we  are  fitting  are  shown  in 
figure  19,  along  with  the  corresponding  experimental  reverse  bias  value. 

Based  on  figure  24b,  we  set  the  Cd  alloy  composition  profile  in  the  top  P+  layer,  x(III),  as  a 
fitting  parameter  in  the  form  of  a  linear  grading  xtop— >0.3  (with  xtop>0.3).  The  two  other  fitting 
parameters  are  the  ones  that  are  found  to  have  the  most  influence  on  the  calculated  device 
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characteristics:  IF1F1I  in  the  absorber  layer,  which  is  related  to  the  magnitude  of  Auger- 1 
processes  and  the  SRH  characteristic  lifetime  rno  =  tpo  (assuming  the  same  for  electrons  and 
holes)  in  the  absorber  layer,  which  is  related  to  the  concentration  of  traps,  impurities, 
dislocations,  or  other  elements  causing  SRH  recombination,  as  well  as  their  capture  cross- 
sections.  We  observe  that  IFiF2I  mostly  has  an  influence  on  the  value  of  (Jraax-Jmin)  and,  therefore, 
the  strength  of  the  NDR.  The  SRH  lifetime  xno  (xpo)  has  an  influence  on  the  values  of  both  Jmm 
and  Jmax,  but  mostly  on  Jmin.  Table  2  summarizes  the  parameters  used  in  our  device  model  to  fit 
the  experimental  I-V  characteristics.  In  this  case,  the  temperature  dependence  of  IFiF2I  and  xno 
(xpo)  is  taken  into  account.  The  Cd  alloy  profile  in  the  top  P+  layer  x(III)  is  constant  with 
temperature.  The  doping  profile  and  device  geometry  are  taken  from  experimental  data  and 
designed  values.  Compositional  grading  between  each  region  interface  is  taken  into  account.  An 
iteration  scheme  described  in  details  elsewhere  is  used  to  obtain  the  best  values  for  each  fitting 
parameter.  The  fitting  for  each  temperature  is  obtained  after  a  maximum  of  2-3  iterations. 

Ultimately,  the  goal  of  this  fitting  is  to  understand  what  mechanisms  are  limiting  the  leakage 
current  just  before  (Jmax)  and  after  (Jmin)  Auger  suppression.  The  total  leakage  current  Jtot  can  be 
expressed  as  the  sum  of  the  electron  diffusion  current  in  the  top  P+  layer,  Jn,diff(P+X  the  hole 
diffusion  current  in  the  bottom  N+  layer,  Jp,diff(N+),  and  the  absorber  layer  current  Jabs: 

J tot  =  J n,diff{p  )+  Jabs  +  ^ p,diff  )■ 


Table  2.  Extracted  fitting  parameters. 


T  (K) 

if,f2i 

^n0— ^p()  (s) 

Jmax  (A/cm  ) 

Jmin  (A/Cm  ) 

Exp. 

Sim. 

Exp. 

Sim. 

120 

0.29 

2.8xl0~7 

0.24 

0.24 

0.22 

0.21 

130 

0.28 

3.3xl0~7 

0.76 

0.75 

0.52 

0.51 

150 

0.26 

3.7xl0~7 

2.60 

2.51 

1.42 

1.49 

200 

0.25 

3.8xl0~7 

4.78 

4.85 

2.49 

2.41 

270 

0.21 

3.5xl0~7 

6.19 

6.12 

3.19 

3.21 

300 

0.20 

3.4xl0~7 

10.96 

10.83 

5.54 

5.52 

As  seen  in  figure  9,  the  absorber  layer  current  is  a  combination  of  diffusion  and  g-r  current, 
depending  on  the  reverse  bias.  However,  it  is  much  more  interesting  to  define  it  as  the  sum  of  an 
Auger- 1  current  in  the  absorber,  JAi,  an  Auger-7  current  in  the  absorber,  JA7,  and  an  SRH  current 
in  the  absorber,  JSrh: 

J abs  ~  J A\  +  J A1  J SRH  • 
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This  way,  we  can  determine  the  recombination  mechanisms  limiting  the  leakage  current  in  the 
HgCdTe  HOT  photodiodes.  Jmin  and  Jmax  are  fitted  at  each  temperature  using  the  set  of  three 
fitting  parameters  described  earlier,  and  we  extract  from  our  simulation  calculations  the 
corresponding  value  for  each  current  component  mentioned  in  these  last  two  equations.  The 
values  for  the  Auger  and  SRH  currents  in  the  absorber  are  determined  by  switching  “on”  and 
“off”  each  recombination  mechanism  separately  in  our  device  model. 

Figure  25a  shows  the  calculated  magnitude  of  Jmax  after  fitting,  along  with  the  various 
magnitudes  of  the  leakage  current  components.  According  to  our  fitting  result,  the  maximum 
current  before  Auger  suppression,  Jmax,  is  limited  by  a  SRH  current  in  the  absorber  at  T=120  K. 
This  SRH  current  obtained  in  our  model  may  be  associated  in  the  experimental  device  with  a 
background  trap  concentration  introduced  during  growth  and/or  processing.  This  SRH  current 
has  a  weak  dependence  with  temperature.  We  point  out  that,  according  to  our  model,  SRH 
recombination  is  the  limiting  mechanism  for  Jmax  up  to  -130-140  K.  At  T=120  K,  both  Auger- 1 
and  Auger-7  currents  in  the  absorber  are  very  low  compared  to  Jsrh  because  the  intrinsic  carrier 
concentration  in  the  absorber  is  still  relatively  low.  However,  both  currents  rapidly  increase  as 
the  temperature  is  increased.  The  Auger-1  current  in  the  absorber  limits  Jmax  above  -130-140  K. 
Finally,  the  hole  diffusion  current  in  the  bottom  N+  layer  stays  relatively  weak  at  all 
temperatures.  The  electron  diffusion  current  in  the  top  P+  layer  has  some  influence  on  Jmax  above 
270  K.  Figure  25b  shows  the  calculated  magnitude  of  Jinm  after  fitting,  along  with  the  different 
magnitudes  of  the  leakage  current  components.  According  to  our  fitting  results,  the  maximum 
current  after  Auger  suppression  (and  before  tunneling  processes),  Jmin,  is  limited  by  a  SRH 
current  in  the  absorber  between  120  K  and  room  temperature.  This  is  expected  in  LWIR  HgCdTe 
and  can  be  reduced  as  the  HgCdTe  growth  and  processing  technology  is  further  optimized.  The 
Auger- 1  current  JAi,  which  was  found  to  limit  Jmax  in  figure  25a,  is  fully  suppressed  as  the 
reverse  bias  is  increased  and  is,  therefore,  not  shown  in  figure  25b.  The  Auger-7  current  in  the 
absorber  JA7  increases  as  the  temperature  is  increased.  Although  it  is  not  found  to  limit  Jmin,  it 
shows  that  the  Auger-7  current  is  only  partially  suppressed  in  the  device.  The  amount  of  Auger- 
7  suppression  is  -70  %  according  to  our  results.  The  hole  diffusion  current  in  the  bottom  N+ 
layer  stays  relatively  weak  at  all  temperatures,  as  observed  when  analyzing  Jmax.  The  electron 
diffusion  current  in  the  top  P+  layer  has  some  influence  on  Jmm  above  270  K.  This  analysis  gives 
us  significant  insight  on  the  limiting  leakage  current  mechanisms  in  experimental  LWIR  HOT 
devices,  both  before  and  after  Auger  suppression.  Dark  current  caused  by  tunneling  mechanisms 
is  expected  to  occur  at  higher  biases  than  the  one  considered  here  and  is,  consequently,  not 
observed  in  our  simulations. 

The  saturation  current  in  our  experimental  device  (Jmin)  is  next  compared  to  that  in  a  standard 
HgCdTe  DLPH  photodiode  simulated  using  our  device  model.  We  also  compare  Jmin  to  the  dark 
current  in  a  simulated  ideal  HOT  photodiode  with  the  same  fixed  cutoff  wavelength  as  the 
experimental  device  (-10.4  pm  at  100  K).  We  plot  the  dark  current  at  200  mV  reverse  bias 
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versus  T  using  our  numerical  model  for  ideal  HOT  and  DLPH  photodiodes  (x(abs)=0.22), 
characterized  by  a  low  n-type  doping  in  the  absorber,  ND(absorber)=  1x10  cm  ,  and  a  much 
longer  xn0=in0=5  ps  than  what  we  observed  experimentally.  As  seen  in  figure  26,  the  fabricated 
HOT  device  gives  a  lower  saturation  current  Jmjn  than  for  the  ideal  DLPH  at  all  temperatures. 
Above  -200  K,  both  ideal  and  experimental  HOT  devices  have  comparable  dark  currents.  Below 
200  K,  the  experimental  HOT  device  has  a  much  larger  saturation  current  than  the  ideal  HOT 
device,  which  shows  the  room  for  improvement  that  can  be  achieved  by  optimizing  the  growth 
and  fabrication  process.  Improvements  are  especially  needed  in  terms  of  lower  doping  levels  in 
the  absorber  layer  and  fewer  traps/dislocations  causing  SRH  recombination.  At  300  K,  factor  of 
-10  lower  current  density  can  be  achieved  with  this  non-equilibrium  mode  of  operation  by  using 
low  p-type  doped  in  the  absorber  layer. 


(a)  (b) 


Figure  25.  (a)  Experimental  values  for  Jmax  versus  temperature.  Also  shown  are  the  calculated  Auger-1  current  in  the 
absorber  JA1,  Auger-7  current  in  the  absorber  JA7,  SRH  current  in  the  absorber  JSRH,  electron  diffusion 
current  in  the  top  P+  layer  Jn,diff(P+).  and  hole  diffusion  current  in  the  bottom  N+  layer  Jp,dtff(N+). 

(b)  Experimental  values  for  Jmin  versus  temperature.  Also  shown  are  the  Auger -7  current  in  the  absorber 
JA7,  SRH  current  in  the  absorber  JSRH,  electron  diffusion  current  in  the  top  P+  layer  Jn  din(P+),  and  hole 
diffusion  current  in  the  bottom  N+  layer  Jp-djff{N+) 
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Figure  26.  Calculated  dark  current  density  at  200  mV  reverse  bias  for  an  ideal  DLPH  and  HOT 
devices  (with  x(abs)=0.22)  and  experimental  Jmin  versus  temperature. 


8.  Conclusion 


We  have  demonstrated  Auger  suppression  devices  at  room  temperature  on  p-type  As  doped  LW- 
HgCdTe  materials  grown  by  MBE.  We  employed  a  three-layer  device  architecture  P'/7t/N+  in 
this  work.  Data  shows  about  a  50%  reduction  in  dark  current  is  achieved  at  300  K.  A  principal 
challenge  was  to  obtain  low  p-type  doping  levels  in  the  ti- region.  This  issue  was  overcome  using 
a  novel,  deep  diffusion  process,  thereby  demonstrating  a  successfully  low-doped  p-type  HgCdTe 
in  MBE  grown  material.  Results  indicate  Auger  suppression  in  P+/7i/N+  devices  above  130  K  and 
the  obtained  saturation  current  is  of  the  order  of  5  A/cm'  on  these  devices  at  300  K.  The 
threshold  voltage  is  450  mV  at  300  K  and  100  mV  at  130  K  for  250-micron  circular  devices. 
Results  show  reducing  the  series  resistance  could  reduce  this  further.  Near-classical  spectral 
responses  were  obtained  at  250  K  and  100  K.  Measured  non-AR  QE  at  100  K  was  0.57.  By 
fitting  the  temperature-dependent  current-voltage  experimental  data  for  one  selected  device 
where  clear  NDR  is  observed,  we  confirmed  that  the  NDR  is  due  to  Auger  suppression.  More 
specifically,  NDR  is  attributed  to  full  suppression  of  Auger- 1  processes  and  partial  suppression 
of  Auger-7  processes.  The  percentage  of  Auger-7  suppression  is  -70  %.  After  Auger 
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suppression,  the  remaining  leakage  current  is  mostly  limited  by  a  high  SRH  recombination 
component  in  these  devices. 

Also  in  this  study,  we  have  demonstrated  a  novel  method  of  shallow  arsenic  ion  implantation  and 
deep  diffusion  anneal  to  obtain  controllable  deep  diffusion  of  low  arsenic  concentrations  in  MBE 
grown  HgCdTe  DLPH  samples,  while  retaining  the  P+/7t  step  structure.  We  achieved  low  arsenic 
concentrations  of  (0.8  to  2)xl0  cm  in  the  absorber  region  using  low  Hg  partial  pressure  and 
long  anneal  times,  Finally,  we  discussed  the  behavior  of  arsenic  in  HgCdTe  as  a  function  of 
variations  in  the  deep  diffusion  parameters, 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


Ag 

silver 

BLIP 

background-limited  performance 

BTB 

band-to-band 

FTIR 

Fourier  transform  infrared 

g'r 

gcncration-=recombi  nation 

HgCdTe 

mercury  cadium  telluride 

HOT 

high-operating-temperature 

InSb 

indium  antimony 

I-V 

current-voltage 

LWIR 

long-wave  infrared 

MBE 

molecular  beam  epitaxy 

NDR 

negative  differential  resistance 

PEC 

process  evaluation  chips 

PV 

photovoltaic 

QE 

quantum  efficiency 

SE 

spectroscopic  ellepsometry 

SIMS 

secondary  ion  mass  spectroscopy 

SRH 

Shockley-Read-Hall 

SWaP 

size,  weight  and  power 

TAT 

trap-assisted  tunneling 
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